A Theory of Inertial Range Similarity in Isotropic Turbulence 
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We consider equilibrium statistics for high Reynolds number isotropic turbulence in an incom- 
pressible flow driven by steady forcing at the largest scale. Motivated by shell model observations, 
we develop a similarity theory for the inertial range from clearly stated assumptions. In the right 
variables, the theory is scaling invariant, but in traditional variables it shows anomalous scaling. 
We obtain the underlying probability density function, the scaling exponents, and the coefficients 
for the structure functions. An inertial range length scale also emerges. 

PACS numbers: 47.27.Eq, 47.52. +j, ll.80.Cr, 47.27.Gs, 05.45.Jn, 47.27.Jv, 47.27.eb, 47.27.ed 
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Turbulence is often called the last great unsolved prob- 
lem in classical physics. Being responsible for rapid mix- 
ing and energy dissipation in fluids, turbulence is also 
critical in many applications. Turbulence research aims 
to describe an ensemble of flow realizations statistically 
using moments or probability density functions (PDFs). 
The simplest and most symmetric problem to consider 
is steady-forced isotropic turbulence in an incompress- 
ible fluid. At high Reynolds numbers, large-scale forc- 
ing is separated from the dissipation by a wide range of 
scales where inertial forces rule and the physics is univer- 
sal. In 1941, Kolmogorov developed the first theory 
for this inertial range. His theory removes scale depen- 
dence by normalization. Unfortunately, his scale invari- 
ance is at odds with observations. The discrepancy is 
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FIG. 1: Compensated exponents. Experimental data: longi- 
tudinal (+) from 0, transverse (x) from Q. Our shell model 
data (o); note the statistical convergence is poor for p > 5. 
The curve is our theory with (3 = 1.83 and C,-2 = —0.82. 
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FIG. 2: Normalized one dimensional PDF obtained from our 
theory, n is the shell number which increases towards smaller 
scales. 



called "anomalous scaling." According to Frisch p. 
ix], "modern work on turbulence focuses to a large ex- 
tent on trying to understand the reasons for the partial 
failure of the 1941 theory." 

In the inertial range, the pth order moment or struc- 
ture function S p for the velocity difference Sv(l) between 
two points is a power law in the separation distance I 
@, i-e., S p cx fa. Kolmogorov 0,0 suggested ( p = p/3. 
This way, S p (l) / S^ 2 (I) is independent of I. Moreover, 
the PDFs for Sv at different / collapse to a single curve 
when plotted in units of standard deviations. For a wide 
inertial range, Kolmogorov's scaling invariance seems an 
obvious expectation. It is, however, refuted by experi- 
mental data |2|, y, |4| . ( p is not a linear function of p (Fig. 
^1 , and the shape of the PDF varies with scale as in Fig. 
121 The deviation of ( p from p/3 is the anomalous scaling. 
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Naturally, the central question is this: Why no scaling 
invariance? We suggest the simple answer that scaling 
invariance emerges in the PDF not for Sv, but rather 
for a different function. This function is easiest to find 
in wave number space. There, k • u = together with 
u(— k) = u*(k) leaves two real variables at each k. 
Equivalently, the complex helical waves decomposition 5] 
expresses u(k) using a right and left-handed amplitude. 
Let £k be their joint PDF in wave number shell k = |k|. 
Using polar form instead of two real amplitudes, we can 
expand azimuthally 

oo 

&(r,0)= ^i m) (r)cos(m(0-# m )). 

m=0 

In shell model statistics, we have found scaling invariance 
of o exp, which is the mean radial profile of on 



a logarithmic abscissa. This paper shows how scaling 
invariance of o exp leads to a theory of anomalous 
scaling. We do not consider other ^ m ' ) here. 

Let us first consider the shell model data. Although 
crude models of the Navier-Stokes equations in wave 
number space, shell models allow us to gather tur- 

bulence statistics at asymptotically high Reynolds num- 
bers. Moreover, shell models exhibit anomalous scaling 
0,0. We use Zimin's shell model 0,E1E3- Built on 
complex helical waves, it inherits the correct balance of 
left- and right-handedness from Navier-Stokes. Wavelets 
partition wave number space into shells 2 n 7r < |k| < 
2 n+1 7r, and, in the spirit of Lorenz |j, two collective 
variables R n and L n model all wavelet coefficients within 
each shell. Using toroidal S n — R n + L n and poloidal 
D n = R n — L n velocities, the evolution equations read 
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where T m = for m < -1, T_ x = 0.1935, T = 0, and 
T m = 2 5 / 2 T_i for m > 1; v is the kinematic viscosity; 
and F = (Fs, n > Fd,u) is the forcing. In polar form, we 
have (S n , D n ) = A„(cosi?„, sini? n ), where A n is the shell 
amplitude. The shell energy is E n — p n A^/4:, and the 
enstrophy Z n = 7r 2 4 n p n A\, where p n = 77r8 ra /18 is the 
wavelet density. With E — and Z — ^Zn, we 

form the Reynolds number R x = (E) ^- 1 (20/(3 (Z))) 1 / 2 . 
Our focus is equilibrium statistics for isotropic turbulence 
in the classical setting of a fixed integral scale (e.g., fixed 
box size) and steady helicity-free forcing at the largest 
scale. Thus, we truncate so n > and set the 

forcing as F s , n — S n o, F D<n = 0. Using VODE 0, we 
generate an equilibrium ensemble of 260,000 realizations 
starting from random large scale initial conditions (Rx — 
3 x 10 6 and vj \F\ = 1.3 x 10" 5 ). 

Let ip n be the PDF of lnA„, i.e. ip n (x)dx = Pr{x < 
lnA„ < x + dx}. We compute ip n from our ensemble 
using kernel density estimation. The result (Fig. |3J) re- 
veals scaling invariance in the inertial range: shifts and 
abscissa-scaling collapse the graphs of ip n (x) / e 2x to a sin- 
gle curve. Thus, one function / generates i\) n by the 
similarity: ip n (x) = e 2x T n f((x-fi n )/a n )/a n . This for- 



mula expresses scaling invariance of £ fe 
be the PDF for A n . Then, cj> n (x)dx - 



(0) 



o exp. Let (f> n 
Vi{x < A n < 



x + dx} = Pr{lnx < ]xiA n < hix + din a;} so that 
4>n{x) — ip n (lnx)/x, and our similarity becomes 



$n(x) _t jl I ( x_ 
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where g = f o ln and jl n = e^™ . 

Now we develop the theory. Motivated by the shell 
model, let us assume (0 and seek those functions g al- 
lowing 

Sp{n) = C p 2-^ n ~ no \ 
We use Mellin transforms, i.e. 



G{z)=m[g{x),z] 



x z 1 g(x)dx, 



to express S p (n) as moments of <f> n . For posi- 
tive constants a and q, operational rules |l3| include 
m[g(ax),z] = a- z G{z) and M[g(x q ), z] = G{z/q)/q. Us- 
ing these, 

S p (n) = P p n /2 (AP) = p p n /2 m[cf >n (x)/x,p + 2] 

= p p n /2 T n pP+ 2 G((p + 2)a n ). (4) 

Using So(n) = 1, we eliminate r n so that S p (n) = 
p p n /2 flP n G((p + 2)<j n )/G(2a n ). Setting p = 3 yields 
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FIG. 3: (Color online) Statistical shell model data. The top 
panel shows the collapse in the inertial range and our theo- 
retical curve (heavy black). 



An = C 3 lA V« 1/2 (G(2a„)/G(5 ( r„)) 1/3 2-ft(»-»o)/3. Con _ 
sequently, 

6p(n) -° 3 {g(5^}) G(2a n ) ■ 

The assumption S p (n) — C p 2~'> !, ( n ~™ ' implies 

( G(2a n ) \ p/3 G((p + 2)a n ) = A{p){n _ no)+B{p) 
\G(5a n )J G(2a n ) 

for some A(p) and B(p). We set no so that B contains 
no multiple of A. With r\ = In o G we obtain 

| (r?(2CT„) - r)(5a n )) + n((p + 2)a n ) - rj (2cr n ) 
= A(p)(n~n ) + B(p). 

Assuming the existence of an invertible function s such 
that s(n) = a n we have w — s(n) n = s^ 1 (w) and a 
linear functional equation for rj: 

| (n(2w) - n(5w)) + n {(p + 2)w) - n (2w) 

= A(p)( S - 1 (w)-n a )+B(p). (5) 

The homogeneous solution is rj(w) — aw + b. It does not 
affect Sp(n). In fact, a = b = means normalization, 
but no loss of generality. The non-homogeneous solution 
is then 7/(11;) = d\w^ + c^lntf, where di, di and (3 are 
constants. Special solutions n = d\ (lnu;) 2 + ^lnw and 



rj = d\w\i\w-\-d2 In w substitute for the degenerate cases 
(3 = and (3 = 1. Shell model data indicate (3 « 1.83, so 
we present analysis for (3 ^ 0, 1. That is, G(z) = e v ^ z ' = 
z d2 e dizl3 . Experimental 01 an d (our) computational ev- 
idence show that the joint PDF is finite and non-zero 
at the origin. Thus, we add the theoretical assumption 
of a finite non-zero value for g(x) = 9Jt~ 1 [z d:! e dlZ ,x] 
as x —* + . This implies c?2 = — 1- Inserting rj = 
d\w^ — lnui into JHJ, we obtain B(p) = — pln(2/5) — 
In(p/2+l), A(p) = di (f (2" - 5") + (p + 2)^ 3 - 2"), and 
<r n = s(n) = {di/di) 1 ^ (n — hq) 1 ^ , where d± is con- 
stant. Expressing d± in terms of Q-2, we find for p > —2, 



C P = (C3 + 2C-2) 



(p/2 + l)^-(p/2 + l) 
(5/2)" -5/2 
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P > -2. 



C-2, (6) 



(7) 



Serving only to normalize a n , the constant d\ is absent 
from C P and C p . We choose di so that (d^/di) 1 ^ = 1/2. 
With k = (C 3 + 3C- 2 /2)/((5/2) /3 - 5/2) we have 



G(z) 



(8) 



All together, we have five parameters: C3, C3j A n o- 
Kolmogorov's 4/5-law implies (3 = 1- To find the oth- 
ers, we use cumulants of In A n . Let v I'„(s) be the char- 
acteristic function of i^ n (x) = e x (f> n (x). By substitut- 
ing x = In it, we obtain ^> n {s) = SDt [<p n (u), is + 1] = 
P~ ls S ls {n) = p- is C is 2-<"( n - n °\ Clearly, ln*„(«) is a 
linear function of n, as are the cumulants (((In A n ) )) = 
l(-id/ds) m ln^ n (s)] s=0 . In particular, 2 ((lnA„)) + 1 = 
-ln(77r/18) - 3nln2 + 21n(G 3 5/2)/3 - {k{(3 - 1) - 
C-a)(n - n ) In 2, and 4«(lnA„) 2 )) - 1 = -n(3{(3 - 
1) (In 2) (n-n ). We plot ((lnA n )) and (((ln^„) 2 )) for 
our shell model data in Fig. 0] The data clearly identifies 
two lines. Using them, we find C3 = 3.85, (-2 = —0.825, 
(3 = 1.83, n = —3.7. These values collapse the data in 
Fig. with 

r n /a n = 2 Pn (5C 3 /2)- 2/3 2-<»— 



In 



(5G 3 /2) 1/3 (2(^- 2 )("- n »)/p„) 1/2 
^-i(n-no) 1 /' 3 . 



As the virtual origin of the inertial range scaling, no iden- 
tifies a turbulence length scale exclusively from inertial 
range properties. 
Our theory yields 
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(|K|ln2)- 1 /' 3 /cr„ 



,(9) 
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FIG. 4: (Color online) Cumulant data used to identify no, (3, 
C3, C-2 for the shell model. 



where numerical techniques |15| provide the inverse 
Mcllin transform. Recall that (f> n (x)/x represents the 
first term in the azimuthal expansion of the joint 
PDF so to compare with one-dimensional PDFs 
from the literature, e.g. 0, 0, 0, we integrate 

(f> n ((S% + D^) 1//2 J over D n . The result is an approxi- 
mation to the PDF for S n (Fig. The shape of the 
PDF differs between shells. It is roughly Gaussian for 
small n — n , but develops a sharp peak and broad tails as 
n increases. These are classical characteristics of anoma- 
lous scaling. 

Theory often requires that £ p increases monotonically 
. This is a consequence of assuming an upper bound for 
the velocity so as to avoid supersonic speeds 0. Corre- 
spondingly, the PDF for the velocity increment can not 
have an infinite tail, but must have compact support. 
Our function g has compact support when (3 < 1, but not 
when (3 > 1. Importantly, our data yields f3 ~ 1.83 > 1. 
Consequently, £ p has a maximum and even decreases to 
negative values when p is large. Negative values of Q, are 
not unphysical. In fact, they occur in 2D turbulence |l6| . 
As a mathematical problem, turbulence described by the 
incompressible Navier-Stokes equations does not know of 
supersonic speeds and has its own statistics. We can not 
impose constraints on that statistics. It is a different 
matter that the incompressible equations incorrectly de- 
scribe fluid flows at supersonic speeds. It is wrong to 



require d(, p /dp > simply because the condition imposes 
compact support on the PDF. In contrast, d 2 ( p /dp 2 < 
follows from ((A 2 n )) > 0. 

Our theory has scaling invariance in the form of the 
similarity formula but also reproduces the well 

known anomalous scaling when expressed in traditional 
variables (Fig. [2J. Thus, by choosing the right vari- 
able, anomalous scaling disappears and we obtain a col- 
lapse across scales as in Figure 3. Our theory rests on 
four inertial range assumptions: scaling invariance of 
£^ o exp (Eqn.®), existence of the exponents £ p , a fi- 
nite £fc°^(0) 7^ 0; an d an invertible a n . These assump- 
tions hold for our shell model. We can only conjecture 
that they also hold for Navier-Stokes. To test our theory, 
one should aim to reproduce the collapse in Fig. [3] That 
is, sample ln(£J(|k|)) for various k — |k| in the inertial 
range, i?(|k|) being the three-dimensional energy. Call 
the PDF ipk and plot the lower part of Fig. with ipk 
replacing i\) n . If curves for different k can be brought to 
coalesce by shifts and horizontal scaling, then our theory 
applies and gives the analytical expression for ipk as well 
as many other formulas. 
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